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Abstract. 

We study cosmological perturbations in the case that present-day matter consists of 
a mixture of inflaton and curvaton decay products. We calculate how the curvaton 
perturbations are transferred to its decay products in the general case when it does 
not behave like dust. Taking into account that the decay products of the inflaton 
can also have perturbations results in an interesting mixture of correlated adiabatic 
and isocurvature perturbations. In particular, negative correlation can improve the fit 
to the CMB data by lowering the angular power in the Sachs- Wolfe plateau without 
changing the peak structure. We do an 11-parameter fit to the WMAP data. We 
find that the best-fit is not the 'concordance model', and that well-fitting models do 
not cluster around the best-fit, so that cosmological parameters cannot be reliably 
estimated. We also find that in our model the mean quadrupole (I — 2) power is 
1(1 + 1)Ci/2tt = 1081^K 2 , much lower than in the pure adiabatic ACDM model, which 
gives 1262/xK 2 . 



PACS numbers: 98.70.Vc, 98.80.-k, 98.80.Cq 
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1. Introduction 

The building blocks of what is today known as the curvaton model were introduced in 
[TJI2] (some early discussion can be found in OH])- The first detailed application to the 
observed CMB anisotropics was done in jH], in the context of the pre-big bang scenario. 
The model was then applied to inflation in detail and given the name "curvaton" in 



In the curvaton model the perturbations of the cosmic microwave background and 
large-scale structure originate from a field, the curvaton, which is different from the 
field(s) driving the dynamics of the very early universe, be it inflation, pre-big bang 
collapse or something else. In what follows, we adopt the language of single field inflation 
for describing the primordial universe, though our considerations could apply to other 
cosmological scenarios as well. 

During inflation and reheating, the curvaton is assumed to be completely 
sub dominant. The mass of the curvaton field is also assumed to be much less than the 
Hubble parameter during inflation, so that the curvaton field value is frozen, and the 
curvaton will acquire a spectrum of quantum fluctuations which become classical. After 
inflation, when the Hubble parameter falls and becomes of the order of the curvaton 
mass, the curvaton will start to oscillate around the minimum of its potential. Assuming 
the potential to be quadratic about the minimum, the energy density in the curvaton 
field, averaged over oscillations, decays like a -3 , where a is the scale factor. Taking the 
universe after reheating to be dominated by radiation, the relative contribution of the 
curvaton to the energy density increases. (Since the curvaton field value and therefore 
energy density is constant before the start of the oscillations, its relative contribution 
increases even faster during that period.) 

We assume that the curvaton decays into baryons, cold dark matter (CDM), leptons 
and photons. The leptons are assumed to consist of (or later decay into, depending on 
when the curvaton decays) electrons and neutrinos. (Electrons play no role in our 
discussion, and will not be mentioned from now on.) If the curvaton decays before it 
contributes to the background, the perturbations inherited by its decay products will be 
of the isocurvature type, whereas if the curvaton completely determines the background 
when it decays, the perturbations will be adiabatic. This conversion of isocurvature 
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perturbations into adiabatic ones is the essence of the curvaton mechanism. If the 
curvaton decays between these two extremes, its decay products will inherit a mixture 
of correlated adiabatic and isocurvature perturbations jHj- This has been studied in 

ia cni en ca us nsi cni en 

In the present paper, we extend in two respects previous studies (building 
particularly on [TH]) of the case that the curvaton decays while it contributes to the 
energy density but does not completely dominate it. 

First, it has usually been assumed that one can treat the curvaton as a pressureless 
fluid in calculating its decay. However, as this description is valid only as an average 
over oscillations, one would not expect it to apply in the case of a rapid decay, or when 
the field decays just as it is starting to oscillate (recently studied in [17J). We will 
calculate the decay without this assumption, using the curvaton field equation, and will 
obtain the limit of validity of the dust approximation. 

Second, it has usually been assumed that even when present-day radiation and 
matter consist of a mixture of inflaton and curvaton decay products, the inflaton decay 
products have no perturbations. We will take into account the possibility that both the 
inflaton and curvaton decay products carry some (uncorrelated) perturbations, resulting 
in an interesting mixture. This case has been studied in EI and is also reminiscent 
of double inflation models (21 EH EH 1201 El 122]- Two limiting cases are the usual inflaton 
scenario, where all present-day matter comes from the inflaton, and the original pre-big 
bang curvaton scenario j3], where everything comes from the curvaton. 

In section El we consider the decay of the curvaton, solve the decay equations 
numerically and discuss recovering the dust approximation. In section El we fit our 
spectrum of perturbations to the measurements made by the WMAP satellite. We 
discuss the interesting features of the spectrum, in particular suppression of the low 
multipoles and negative running of the spectral index. Having two independent sources 
of perturbations allows for a spectrum where one can have different behaviour of the 
low and high multipoles. The acoustic peak structure can remain very close to the 
pure adiabatic case, while the angular power on the Sachs- Wolfe plateau can be much 
reduced, leading to better overall fit than the pure adiabatic ACDM model. We also 
find that the inclusion of isocurvature modes opens up the space of well-fitting models, 
so that one cannot determine the cosmological parameters from the WMAP data alone. 
A detailed study of the estimation of the isocurvature contribution and the cosmological 
parameters with our perturbation spectrum will be presented later 

2. The curvaton decay 

2.1. The set-up 

We assume that there is a bath of baryons, CDM, neutrinos and photons produced 
by the decay of the inflaton field. We take the perturbations in these fluids to be 
adiabatic, as expected from single field models of inflation j21J. We also assume that 
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there is a curvaton field which is sub dominant, but frozen at some field value so that its 
relative contribution to the energy density rises. The curvaton field is assumed to have 
a spectrum of perturbations, acquired during inflation, which is uncorrelated with that 
of the inflaton decay products. 

Like the inflaton, the curvaton is assumed to decay into baryons, CDM, neutrinos 
and photons, which inherit the perturbations of the curvaton. The final state will 
therefore consist of a mixture of products of inflaton and curvaton decay, both with 
their own spectrum of perturbations. 

We will follow the system numerically from the time that the curvaton is frozen and 
subdominant until the time it has completely decayed. We will see how the resulting 
spectrum of perturbations in the decay products depends on the initial conditions and 
parameters of the model. Our analysis of the decay is a generalisation of that in |16j . 
where the curvaton decay was followed numerically assuming that the curvaton behaves 

pressureless fluid. 



2.2. The equations 

The background. The spacetime is a perturbed Friedmann- Robert son- Walker universe. 
In accordance with the generic inflationary prediction, we take the background spacetime 
to be spatially flat. In the uniform curvature gauge, the metric reads (we consider only 
scalar perturbations) 

ds 2 = -(1 + 20)dt 2 + 2aB, l dtdx i + a 2 S ij dx i dx j . (1) 

Throughout, we follow the notational conventions of ^21E3I2I3; except that we 
prefer to use TZ, the curvature perturbation on the comoving hypersurface, rather 
than (, the curvature perturbation on the uniform energy density hypersurface. On 
superhorizon scales, the two are related by 1Z = 

As sources, we have radiation (r), matter (m)\ and the curvaton field (a). The 
relation between the background sources and geometry is given by 

p + 3H(p + P) = , (3) 

where H = a /a, Mpi = 1/ ^8irG-^ is the (reduced) Planck mass and p and P are the 
total energy density and pressure, respectively. In terms of the individual components, 
we have 

^aPa = P , ^aPa = P , (4) 

where a = r,m, a. For radiation and matter we have P r = |p r , P m = 0, so that if there 
was no energy transfer between the components, their energy densities would decay like 

| The identification of radiation and matter as baryons, CDM, neutrinos and photons will be discussed 
in section mi 
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p r oc a 4 and p m oc a 3 . Since the curvaton decay transfers energy into radiation and 
matter, the behaviour of the energy densities is given instead by 

p a + 3H( Pa + P a ) = Q a , (5) 

where Q a is the energy density transfer per unit time to the fluid a. 

We take the curvaton field to be minimally coupled to gravity and to have the 
potential V(a) = \m 2 a 2 . We model the curvaton decay phenomenologically by 
introducing constant decay rates to both radiation and matter into the equation of 
motion, as in ^H]- The curvaton equation of motion for the background then reads 

a + (3H + T)a + m 2 a = , (6) 

where a is the background value of the curvaton field and r = T r + T m is the curvaton 
decay rate, with T r and T m being constants. Using this effective description of decay, 
we ignore the possibility of parametric resonance j2Hl I2Z| , as is usually done for curvaton 
models. Thermal corrections to the potential and to the decay rate [2H] could also be 
important, since the curvaton is coupled to the thermal bath of matter and radiation 
[2H]; we ignore these too. 

The energy density and pressure of the curvaton are given by 

■ 2 i 2 2 

P*=2 a + 2 m a 

P a = -a 2 - -m 2 a 2 , (7) 
and the energy transfers are 

Q r = T r (p a + P a ) 

Qm ^mi^Pcr P<r) 

Q a = -(T r + T m )( P(T + P a ) . (8) 

Note that p a + P a = a 2 . In [THj, it was noted that we have an upper limit T rn /T r < 10 -6 
from the requirement that the curvaton decay should be complete before big bang 
nucleosynthesis (BBN) at z ~ 10 10 and matter should be subdominant until the matter- 
radiation equality at z < 10 4 . This limit only holds if all matter and radiation come 
from the curvaton, which is not true in our case (or in the case of [USD- We will discuss 
the limit on T m /T r in section 

For comparison with the dust treatment of the curvaton in JH], let us define |30j : 

Pa + P a = (7 + lp) Pa ) ( 9 ) 

where 7p CT and ^ v p a are the average and the periodic part, respectively, of p a + P a over 
an oscillation. In the regime where the curvaton oscillation is much more rapid than the 
decay, m ^> T, we can neglect j p . Moreover, for a quadratic potential 7 = 1, so that in 
the rapid oscillation limit we approximately recover the expressions of [To] : 

Qr f"VPcr 

Qrn = FmPa ■ (10) 

However, we do not expect (|10|) to be a good approximation when the field does not 
oscillate rapidly when it decays. 
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The perturbations. We are interested in the behaviour of perturbations only in the 
large-scale limit, when spatial gradients can be neglected. Then, the equations for the 
perturbations read (after eliminating the metric perturbation 0) ^21 QBj 

5p a + 3H(5p a + 5P a ) = 5Q a - ^Q a , (11) 

zp 

where 5p a ,5P a and 5Q a are, respectively, the perturbation in the energy density and 
pressure of, and energy transfer to, the component a. The total perturbation obeys the 
covariant conservation equation 

5p + 3H(5p + 5P) = , (12) 

where 

5p = T, a 5p a , 5P = Y* a 5P a . (13) 
The equation of motion for the curvaton field perturbation 5a is 



5a + (3H + T)5a + m 2 5a = (2m 2 a + Ta) ^ + 3Ha - -) ^ . 
The perturbations in the energy density and pressure of the curvaton are given by 



(14) 



t- • r • 2 r ■ 2 

5p a = a5a + m ada H a 

2p 

5P a = &5a - m 2 a5a + —a 2 , (15) 

2p 

and the perturbed energy transfers are 

5Q r = T r {5 Pa + 5P a ) 
5Q m = T m (5p a + 5P a ) 

5Q a = -(T r + T m )(5p a + 5P a ) . (16) 
We have assumed that the decay rates are constant also at the perturbed level, as in 



2. 3. Behaviour of the background 

The parameters of the background are the curvaton mass m, the decay rate T, the 
ratio r m /r r , as well as the initial values which consist of the initial curvaton field value 
o"o and the initial radiation energy density p r Q. Under the assumption that matter is 
subdominant, the initial energy density of matter makes no difference. Since we assume 
that initially H ^> m,T, we see from (jnj) that the curvaton is essentially frozen at 
a constant value to which it has been driven during inflation, so its initial velocity 
is negligible (though non-zero). For the curvaton to be subdominant, we must have 
p^ > \m 2 al. 

When the Hubble parameter becomes of the order of the curvaton mass, H ~ m, the 
curvaton rolls down the potential and starts to oscillate. When the Hubble parameter 
becomes of the order of the decay rate, H ~ T, the curvaton decays into radiation and 
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matter. However, since the energy transfer in (JSJ) is proportional to a 2 , there is no decay 
as long as the field is frozen. This means that if m < T, the decay will not start at 
H ~ r, but only later at H ~ m, when the field will simultaneously roll down the 
potential and decay. So, the condition for curvaton decay is H ~ min(m, V). Note that 
this makes it a natural possibility (in the context of our phenomenological treatment) 
for the curvaton to decay already when it starts to oscillate. 

If cr ^ Mpi, we see from © and © that the curvaton starts contributing 
significantly to the energy density before it starts to oscillate. For o"o 3> Mpi, a period 
of curvaton-driven inflation will ensue and the inflaton-curvaton model will in fact be a 
double inflation model. In a realistic model, the curvaton potential is not expected to 
be exactly quadratic due to Planck scale -suppressed non-renormalisable terms, so we 
should not apply our calculation to the case <To 3> Mpi. At any rate, we want to study the 
case when both the inflaton and curvaton decay products contribute to the observed 
radiation and matter, so the region of parameter space where there is a long second 
period of inflation is of no interest to us. Note that if in Planck-scale inflation the values 
of fields such as the curvaton randomly sample a more-or-less uniform distribution, a 
value o"o ~ Mpi is much more likely than a value o~q <C Mp\. For values o"o ~ Mpi and 
T > m, the curvaton will decay as it is starting to contribute to the energy density, 
making it a not unnatural possibility that both the inflaton and curvaton contribute to 
the radiation and matter density. 



2.4- Behaviour of the perturbations 

In the beginning, the system consists of radiation and matter from the inflaton 
decay, carrying one spectrum of perturbations, and the curvaton, carrying another, 
uncorrelated, spectrum of perturbations. In the final state after curvaton decay, only 
radiation and matter, carrying a mixture of the two spectra, are left. 

The curvature perturbation of the component a is (in the uniform curvature gauge) 

n a = h 5 ^ , (17) 

Pa 

and the total curvature perturbation is 

TZr — E a --— 7?. Q = H-— . (18) 
P P 

We are interested in how the perturbations of the curvaton are transferred to its 
decay products. Usually, the curvaton is assumed to behave like dust, P a = 0. Since the 
equation of state is then (trivially) barotropic, P a = P a (p a ), the curvature perturbation 
TZ a is conserved on superhorizon scales [23- Then the ratio TZ^/IZ™ (using only the 
part of IZ™ 1 that comes from the curvaton decay) gives a useful measure of how the 
curvaton perturbations are transferred to the fluid a. When the curvaton is treated like 
dust we have, from © and |T7|) . 
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However, when we treat the curvaton as a field and not as dust, it does not have a 
barotropic equation of state, as is clear from (JJJ), and as a result lZ a is not conserved. 
In this case the curvaton perturbation 7Z a is, from ©, 0, © and (JEj), 

n ° - - H wnrw ■ < 20 > 

Before the curvaton starts to roll down the potential, the field is nearly frozen 
and a is small. As the curvaton starts to oscillate, the absolute value of a increases 
significantly, which translates into a large decrease of lZ a . (When the curvaton field 
oscillates, 1Z a oscillates around a fixed value, so treating it as constant is then justified, 
when averaged over oscillations.) 

Therefore, giving the initial curvaton perturbations in terms of 1Z™ and calculating 
the ratios IZ^/IZ™ is unenlightening, because the results depend strongly on when 
the initial conditions are given. To meaningfully measure the transfer of the curvaton 
perturbations to its decay products, it is better to consider some conserved quantity. 
We therefore define a convenient perturbation variable £ CT by saying that it has the same 
form (fTT?|) as in the dust case, 

k = ~ H TuT^v\ — ' ( 21 ) 
(3H + Y)p a 

but with p a and 5p a given by the expressions for the field case, (JZJ) and (j!5jl . This 
quantity is constant when the curvaton field is frozen and H ^> Y (as we assume initially 
to be the case). The initial value of £ CT is essentially the initial value of Sp a / (3p a ), so 
it is a useful measure of the curvaton perturbations. We will denote the initial value of 
£ by £™ = 72.2, and will be interested in the ratios 72-° w */72.2- The quantity 1Z 2 is related 
to the initial curvaton curvature perturbation by 

K, a - ~2 1^2 
°o 

1 / m 2 an \ „ , v 

= 2( 1 + ^f) R - < 22 » 

where the subscript refers to the initial value. We have traded the strong dependence 
on a 2 for a large normalisation factor of the initial perturbation. 
So, in the beginning we have the following perturbations 

7C = ^1 

C = ^2 , (23) 

where is the spectrum of perturbations inherited by the inflaton decay products, 
which is uncorrelated with In the beginning, the universe is radiation-dominated, 
so the total curvature perturbation is lZ m = 1Z™ = TZ\. 
After the end of curvaton decay we have 

n° r ut = (1 - Ar)^ + A r x r n 2 (24) 

K™ 1 = (1 - A m )n x + A m \ m K 2 . (25) 
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As noted earlier, the curvaton decay has to be complete before the onset of BBN, so 
the universe is still radiation-dominated. The total curvature perturbation and the 
matter-radiation isocurvature perturbation S™1 = — 3(7Z m — 1Z r ) are then 

n out = ^ _ A ^ Ui + ArXr <Ji 2 (26) 
S% = 3(A m - Ar)K x - 3(A m A m , - A r \ r )1Z 2 . (27) 

In (J24J) . ()25|). (|26j) and (|27j). there are two different kinds of factors relating the final 
perturbations to the initial ones. First, the dilution factors A a give the percentage of 
the fluid a coming from the curvaton, 

A a = , (28) 

Pal + Pa2 

where p a \ and p a 2 are the energy densities of the parts of the fluid a that come from 
the inflaton and the curvaton decay, respectively, evaluated in the final state after the 
end of curvaton decay. The value A a = corresponds to fluid a being produced only in 
the inflaton decay and the other extreme A a = 1 corresponds to a being produced only 
in the curvaton decay. 

Second, the coefficients X a tell how efficiently the curvaton decay transfers the 
original curvaton perturbation to the decay products. (Similar coefficients exist for 
the inflaton decay products, but they are identical for radiation and matter since the 
perturbations are adiabatic, and have been absorbed into TZ\.) If everything came from 
the curvaton, A r = A m = 1, we would have X a = IZ^/C". I n other words, multiplying 
by A a translates the initial curvaton perturbations £* n into the final radiation and matter 
perturbations, and the ratio A r / A m measures the relative efficiency of the transfer. 

In section l2~2l we noted that when all radiation and matter come from the curvaton, 
the decay rates have to obey the limit T m /T r < 10~ 6 . In the general case, this relation 
is modified by the dilution factors as follows: 
10 - 6 > 

~ p r 

A-r Pm2 
A m Pr2 

- X~ , (29) 

■rim J- r 

where the energy densities are evaluated at the end of curvaton decay. For A m = A r = 1 
we obtain the limit r m /T r < 10 -6 . Decreasing the matter dilution factor A m only makes 
the limit more stringent, but decreasing A r relaxes the limit, so that T m /r r can be 
arbitrarily large. 

In curvaton models the curvaton spectrum can be non-Gaussian CHI HH UH] • 
This is because usually TZ\ = in (J2IJ), so that if A r is small, \ r lZ2 needs to be large 
and the square of the perturbation can become important. In our model with IZi ^ 0, 
non-Gaussianity would appear in the analogous case where we would insist that the 
inflaton or curvaton perturbations make a sizeable contribution to lZ° ut even though A r 
is close to 1 or 0, respectively. (Likewise for TV^f and A m .) Instead, we assume that the 
perturbations are always small so that non-Gaussianity is negligible. 
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The case of [HI] is recovered from ([24)1. (|25|). (|26|) and (|27|) when there is no matter 
from the inflaton decay (A m = 1), the inflaton decay products have no perturbations 
(TZi = 0) and matter inherits the curvaton perturbations with an efficiency of one 
(A m = 1). Then the only remaining parameter is, in the notation of |16j . r = A r X r . 
The object of the curvaton decay calculation to be presented is to get the range of 
the coefficients A r and A. m in our more general case, and find their dependence on the 
parameters of the model. 

2.5. The decay calculation 

The relevant parameters. We will solve the system of background equations (j2J), JSJ), 
©, © and ((HI) and perturbation equations (JUJ), (JHJ), (fTKJ) and ffTTTf) numerically, using 
the conservation equations © and ()12j) as a check on the calculation. 

Since the calculation is linear in the perturbations, the density perturbations that 
radiation and matter have inherited from the inflaton evolve independently from those 
that they inherit from the curvaton. In calculating how the curvaton perturbations are 
transferred to radiation and matter, we can therefore put the inflaton perturbations to 
zero. 

Since the initial velocity of the curvaton field perturbation is negligible, and the 
initial field perturbation itself only sets the amplitude of the perturbations, the result 
of the perturbation calculation depends only on the background dynamics, as in |16j . 
(This is quite generally the case, because the perturbations are governed by linear second 
order equations which contain no new parameters apart from the initial conditions, and 
the initial velocity is set to zero.) 

Matter is initially subdominant, so we can put the background energy density of 
the matter originating from the inflaton (i.e. the initial matter energy density) to zero. 
The initial velocity of the curvaton field is negligible, so the result of the perturbation 
calculation depends only on the parameters m, T, T m /T r , <r and p r0 . Since the field is 
stuck until H ~ m, the initial energy density of radiation makes no difference, as long 
as it dominates. A larger p r0 simply means that we start at an earlier time, and have 
to wait longer for the curvaton field to become active. 

So, the perturbation calculation depends only on the parameters m, T, T m /T r and 
o"cb specifically on their dimensionless combinations. In addition to the ratio T m /T r , 
we take a / \\f?>Mpi) and T/m as our two independent parameters to vary. In [THj . 
the dynamics boiled down to a single parameter in addition to T m /T r (note that T 
was normalised to the initial value of H, so that a given numerical value of T did 
not correspond to a given theoretical model, but depended on the initial conditions). 
The parameter o"o/(a/3Mp;) controls whether the curvaton dominates before starting 
to oscillate (the case for a / (^/3M Pl ) > 1) or not, and T/m determines whether the 
curvaton decays already as it is starting to oscillate (the case for T/m > 1) or later 
when oscillating. 

In the region <t / (V^Mpi) <S l,r/m < 1, the curvaton oscillates for a long time 
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before dominating or decaying, and the oscillation is rapid compared to the decay. The 
dust approximation should therefore be valid, and we expect to recover the results of 
[To] . In contrast, for Y/m > 1, the curvaton has no time to oscillate before decaying, 
so we expect the results to be different from the dust case (in terms of Q, 7 P is not 
negligible). 

Comparison to the dust case. From the matter and the curvaton, one can form a 
composite quantity which is covariantly conserved for the background [To] : 

r 

m 

Pcomp Pm ~r -p Pa 

p - cm) 

1 comp Y \ / 

When the curvaton behaves like dust, this quantity has a (trivially) barotropic 
equation of state, P CO mp = Pcomp(Pcomp) = 0, so the corresponding curvature perturbation 
is conserved. To study the goodness of the dust approximation, we therefore introduce 
the variable 

t dPcomp 

C,comp — „ ; l^ 1 / 

" Pcomp 

using in analogy with (|2T|) the letter £ to refer to a quantity which is a conserved 
curvature perturbation only when the curvaton behaves like dust. Before the curvaton 
decays, we have p m = (recall that in this calculation we have no matter from the 
inflaton), so that £ comp = £<r, and after the decay is over we have p a = 0, so that 
icomp = T^m- So, i com p is conserved both before and after the decay. The ratio of the 
final and initial values is £™Lp/ £lomp = T^mVC™ = ^m- As noted, when the curvaton 
behaves like dust, £ CT is conserved throughout and we have A m = 1 as found in [To] . 
The value of X m measures the conservation of ^ comp and thus the goodness of the dust 
approximation. 

The numerical calculation. We vary o"o/(v3Mp/) in the range [10~ 4 , 10], going from the 
curvaton being completely subdominant to completely dominating when it is starting 
to oscillate (but avoid a long period of curvaton-driven inflation). We vary Y/m in 
the range [10 -6 , 10], going from the curvaton decaying deep in the oscillating regime to 
decaying just as it is starting to oscillate. We take 100 logarithmically spaced steps for 
each parameter, covering 10 000 models in all. We start our integration when H/m ^> 1, 
the universe is dominated by radiation and the curvaton is frozen at the initial field value 
o"o + 5o~q. For a given value of the parameters we evolve the system until the curvaton 
has decayed, and we are left with radiation and matter only. 

We have studied a number of the above models with different values of the ratio 
r m /r r in the range [10~ 6 , 10 6 ]. The results depend only weakly on the ratio, so we have 
chosen not to do a systematic scan like for the other two parameters. The behaviour 
shown in the pictures below would be qualitatively the same for any value of r m /r r in 
the above range. 
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Figure 1. Background energy densities (up) and £ CO m P (down) for a large initial value 
of the field. On the left, a model where the field decays just as it is starting to oscillate 
and on the right, a model where the field oscillates before decaying. The dotted vertical 
lines mark the nominal start of oscillation and decay. 



There are two reasons why the ratio T m /T r doesn't make much difference to the 
perturbations. First, the behaviour of the background is not affected by how much 
matter the decay produces since we assume that matter is completely subdominant 
throughout. Second, changing T m /T r changes the relative production of radiation and 
matter for both the background energy densities p a and the perturbations Sp a , so the 
perturbation variables which are given by their ratios (for example, H™* = —Sp m / (3p m )) 
are only weakly affected. 

The results of the calculation are shown in Figures to EJ In Figures and [2] we 
show the evolution of the background energy densities and the perturbation variable 
£comp- in Figure H we have chosen a large initial field value, so the curvaton contributes 
significantly to the energy density when it decays. On the left, the curvaton has a 
large decay rate, so that it decays before really starting to oscillate. On the right, the 
decay rate is small, so that the curvaton oscillates rapidly at the time of decay and thus 
resembles a pressureless fluid more closely. In Figure EJ we show the analogous plots 
for a small initial field value, when the curvaton does not contribute significantly to the 
background before starting to oscillate. Note that even in the model on the left where 
the curvaton decays before properly oscillating, we recover the dust result A m = 1. 

The surprising result that the dust approximation is valid even when the field does 



Correlated isocurvature perturbations from mixed inflaton-curvaton decay 
o- /(a/3M P i) = KT 2 , r/m = 1 o- /(V3Mpi) = 10 2 , T/m = 10 7 



13 



0.1 



10' 



10 7 



0.1 



10' 



10 3 10 5 

H/rn 



10 7 




Figure 2. Background energy densities (up) and £ CO mp (down) for a small initial value 
of the field. 

not oscillate when it decays is a general feature, as seen from Figure El which shows 
X m for the whole of our parameter space. The value of A m is practically independent 
of the decay rate T/m: as long as the initial field value a®/ (\f3Mpi) is less than 0.1 
(so that the curvaton does not dominate before oscillating), A m is very close to 1, even 
when the field decays without any oscillations. On the other hand, for larger values 
of a /(\^3Mpi), A m rises rapidly, up to a maximum of about 350 at the edge of our 
parameter space. 

For fitting our spectrum (}26j) . (|2?j) to the data, we need to know, in addition to A m , 
the ratio A r / A m which measures the relative efficiency of the transfer of the curvaton 
perturbations to radiation and matter. In Figure 0] we show A r /A m for our parameter 
space. Like A m , the ratio A r /A m is more sensitive to the initial field value ao/(\^3Mpi) 
than to the decay rate T/m, though it does grow with decreasing T/m. For small values 
of the field, the ratio approaches rapidly, and for large values it tends to 1. The upper 
limit is easy to understand: for large do, the curvaton completely dominates when it 
decays, so its perturbations are adiabatic, and the perturbations of its decay products 
have to be adiabatic as well [23], so radiation and matter inherit the same perturbations. 

We have studied the validity of the dust approximation, and have surprisingly found 
that its results for the decay are valid even when the field does not oscillate when it 
decays. We have also found the range of the parameters A m and A r / A m , which we need 
for fitting our spectrum to CMB data. 
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Figure 3. The matter decay coefficient A m on a logarithmic scale. 




Figure 4. The ratio of the decay coefficients for radiation and matter, A r /A 
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3. The fit to WMAP data 

3.1. The spectrum of perturbations 

Having found the ranges for the parameters for the mixture of inflaton and curvaton 
perturbations in the spectrum (j2ZJ), let us now fit the spectrum to the WMAP 

data and see what such a mixture implies for the CMB anisotropies. 

We have not yet specified the inflaton and curvaton spectra TZi and 1Z 2 (in this 
section we introduce a hat to denote random variable). As the perturbations are assumed 
to be generated during inflation, the spectra depend on the inflationary model. We will 
not choose a particular model, but assume that the spectra are given by simple power 
laws, 



i,) 2 n 2 = N 2 ~k^a 2 (k) , (32) 

where N\ and N 2 are constant amplitudes, n± and n 2 are constant spectral indices and 
k = k/ko with ko = 0.05 Mpc -1 is the dimensionless wavenumber. We have denoted by 
a\(k) and d 2 (k) Gaussian random variables which obey 

(a!> = (a 2 > = , (a t (k)a*(k')) = 6yd® (k - k') . (33) 

If the perturbation spectra were generated during a period of slow-roll inflation, 
the spectral indices would be (to first order in the slow-roll parameters) (HI E2] 

nx — 1 + 2r]i — 6en 

n 2 = l + 2 V2 - 2e H , (34) 

where en = —H/H 2 , and rji and r] 2 are the slow-roll parameters associated with the 
curvature of the potential in the inflaton and curvaton directions, respectively, 

1 d 2 V 



m = m 2 pi 



V 



2 



o 1 d 2 V m 2 

* = M -' v = IIP • (35) 

where (f is the inflaton field and V is the inflaton-curvaton-potential V(o~, ip) = 
^m 2 a 2 +V(ip) (we neglect inflaton-curvaton interactions). Since the mass of the curvaton 
is assumed to be much smaller than the Hubble parameter during inflation, % C 1, and 
since during slow-roll typically \H\ H 2 , one might expect n 2 to be only marginally 
smaller than 1, though ri\ could be more different from (and possibly larger than) 1 due 
to a non-negligible r\\. 

However, in the curvaton scenario there is no need to assume that inflation is slow- 
roll or even that inflationary dynamics are driven by a scalar field [331133 EH]- Also, 
even in single-field inflation the spectral index can be very different from one, while 
having little running [371 EB1 EH] • We do not limit ourselves to slow-roll inflation, and 
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keep the spectral indices as free parameters with a wide range (and do not introduce 
running for the indices). 

In section we calculated the decay of the curvaton into two components, radiation 
and matter, which inherited the curvaton spectrum with different efficiency coefficients 
A Q . In the fit to the CMB, we have four components: baryons, CDM, (massive) neutrinos 
and photons. 

The CDM is assumed to be non-relativistic during the decay, and will therefore 
inherit the curvaton spectrum with the efficiency A m , so its spectrum of perturbations 
is given by (|25|). The photons will inherit the curvaton perturbations with the efficiency 
A r , and have the spectrum (|24jl. During curvaton decay, neutrinos are also relativistic, 
but since the neutrino dilution factor A v can be different from the photon dilution factor 
Ay, neutrinos and photons could get a different spectrum of perturbations. However, 
we assume that there is no net lepton number, in which case the neutrinos and photons 
end up having the same spectrum by virtue of being in thermal equilibrium [TTH |4"U1 |4"T] . 
so there are no neutrino isocurvature perturbations. 

Whether the baryons inherit the curvaton spectrum with the efficiency factor for 
radiation or matter depends on whether the degrees of freedom carrying baryon number 
are relativistic or non-relativistic during curvaton decay. In any case, there is no reason 
for the baryon dilution factor A b to be the same as the CDM or photon dilution factors, 
so there will in general be baryon-CDM and baryon-photon isocurvature perturbations. 
However, since baryon and CDM isocurvature perturbations behave similarly [TBI !4Uj . 
the isocurvature perturbation between baryons and CDM can be transformed to zero, 
at the price of adding more isocurvature between CDM and photons. Then the only 
component with isocurvature perturbations is the CDM, and the spectrum reads 

(<&) ~ 2 K Ut = C 1 - A^Nj^a^k) + A^Njk^a^k) 
(&) ~ 2 S%l a = 3[(4*n - A) + x b(A - A^N^a^k) 

~ ri2 — 1 

- 3[(A cdm X m - AyX r ) + x b (A b X b - A 1 \ r )]N 2 k— S 2 (fc) , (36) 

where x b = p b /p c dm and A& is the baryon decay coefficient, which can be either A r or A m . 

The most general form of a spectrum with an adiabatic plus one isocurvature mode, 
both of which are a combination of two power laws, has six independent parameters: two 
spectral indices and four amplitudes (one for each of the two different powers of k for the 
adiabatic and the isocurvature mode). The expression (|3*6*j) is therefore overdetermined, 
and though the theoretical interpretation of the parameters is clear, their relation to 
observable features is murky. We will write the spectrum ()36|) in a form more suited for 

comparison with observations: 

i 

{^Y^T = N n (cos fk^a^k) + sm0fc^a 2 (fc)) 

(&yS°cdL,-y = N n f s (cos 9 k^a^k)- sin 9 ~k^a 2 (k)^ , (37) 
where the new constant parameters Nn, fs, 4>, & can be written in terms of the old 
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parameters in (J36|) (the inverse is obviously not true, since (|36|) has more parameters 
than flUTj)). 

The new parameters have a transparent interpretation in terms of observable 
quantities: N-r is the amplitude of the adiabatic perturbations, the isocurvature fraction 
fs measures the amplitude of the isocurvature perturbations relative to the adiabatic 
ones, and the angles <fi and 9 determine the proportion of the adiabatic and isocurvature 
perturbations, respectively, that have the spectral index n\ or n 2 - 

The generality of the spectrum (jHZJ) goes beyond the inflaton-curvaton model that 
we have discussed. It is the most general form for the case when baryons, CDM and 
photons inherit power-law perturbations from two independent sources. For example, 
it could describe a double inflation model§. Also, in realistic supersymmetric GUT 
theories, cosmic strings are typically produced after inflation, leading to isocurvature 
perturbations jlHj. The isocurvature perturbations produced by cosmic strings are 
usually considered to be uncorrelated with the pre-existing adiabatic perturbations from 
the inflationary era, but if they were correlated and the strings would decay, a spectrum 
like (}3"T|) might result. 

In the general case, the ranges for the parameters in (jHZj) are as follows. The 
amplitude is positive, the isocurvature fraction fs can be negative, positive or zero, 
the angle <fi is from the interval [0, 7r/2] and the angle 9 is from the interval [0, 7r]. 

Note that the roles of the inflaton and the curvaton are completely symmetric: 
the spectrum ()37|) is invariant under the transformation n\ n 2 , cos<p sin</>, 
cos# sign(cos#) sin# and fs — > — sign (cos ^/g. To avoid scanning the parameter 
space twice, we break the symmetry by writing ri\ = n 2 + An and studying only models 
where the inflaton spectral index is larger than or equal to the curvaton spectral index, 
An > 0. 

If the transfer of curvaton perturbations to matter and radiation were equally 
efficient, A r = X m , then sin9 and cos# would necessarily have the same sign, so that 9 
would be restricted to the range [0,7r/2]. In our model the transfer to matter is more 
efficient than the transfer to radiation, A r < X m , so that 9 has the full range [0, ir] 
when fs > 0, but is restricted to the range [0, 7r/2] when f s < 0. So, the role of the 
decay coefficients \ a is to expand the range of 9 in the spectrum (pTfjl . In our analysis 
we will allow 9 to vary between [0, ir], and at the end remove the forbidden region 
fs<0,9>7r/2. 

§ Note that even in the slow-roll case l|34|) the spectrum l|37ll does not satisfy the consistency relations 
presented in f° r two-field inflation, except when cos 9 = 0. In that case, either the inflaton 
perturbations are zero, or all fluids come from the inflaton and curvaton with the same ratio A a 
(barring a fortuitous cancellation between baryons and CDM in (J2HJ)). 
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3.2. Angular power and (anti) correlation 

The angular power spectra. The temperature-temperature (TT) angular power 
spectrum resulting from the primordial spectra (jHZJ) is 

,3 



cr=^ J f^Mn(k)m+g?A k w k )\ 2 } 

+g^(k)gl s (k)(1z(k)s(kr + n(kys(k))) 

= 4vr J [gf >n (k) 2 (cos 2 k^ 1 + sin 2 P 2 " 1 ) 

+ fl gfA k ) 2 ( cos2 6 ^ ni_1 + sin2 ^ n2_1 ) 

+2gl n (k)f s gf tS (k) (cos cos A;" 1 ^ 1 - sin sin fc" 2 " 1 ) (3? 



dfc 



4vr / ^iV 2 



711—1 



fc 

+ (^(A ; )sin0-/ <s ^(A ; )sin^) 2 r 2 - 1 ] , (39) 

where gfn(k) and gfg(k) are the functions which transfer the primordial curvature and 
entropy perturbations of the mode with the wavenumber k from the early radiation- 
dominated era to the present temperature fluctuation at the multipole I. Written like 
this, the role of fs in scaling the amplitude of the isocurvature modes via gf s (k) is 
particularly transparent. 

For the temperature-polarisation (TE) cross-correlation power spectrum we have 

/Ah £.3 
y 2^<fe(W fe ) +9ls{k)S{k)Y x [g^(k)TZ(k) + g^ s (k)S(k)}) , (40) 

where gfn(k) and gf s {k) are the transfer functions for the polarisation i?-mode (for 
details see e.g. |4*4*]). 

There are six different terms in (J38|) . The first two are due to the primordial 
adiabatic perturbations TZ° ut , and they will be denoted by Cf dl1 and Cf dl2 . The next 
two terms are due to the primordial isocurvature perturbations , and they will 

be denoted by Cf° l and Cf° 2 . Finally, the last two terms arise from the correlation 
between adiabatic and isocurvature perturbations, and will be denoted by Cf orl and 
Cf or2 . The same notation will be used for the components of the TE power spectrum. 

The role of CDM isocurvature perturbations. The acoustic peak structure in the 
measured angular power spectrum looks distinctly adiabatic. CDM isocurvature 
perturbations have valleys where adiabatic perturbations have peaks and are also 
damped faster with increasing I. Therefore pure CDM isocurvature perturbations cannot 
fit the data and are ruled out [15! • However, a mixture where adiabatic perturbations 
dominate the peak structure and CDM isocurvature perturbations give a significant 
contribution at low multipoles remains an interesting possibility. 
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At low multipoles the Sachs- Wolfe effect gives gf s {k) = —2(1 — f v )gfn(k)\\. The 
rapid damping of gf s (k) with increasing I means that if the primordial adiabatic 
and isocurvature perturbation amplitudes are of the same order, the isocurvature 
perturbations modify only the low-/ part (practically the Sachs- Wolfe plateau) of the 
final Ci spectrum observed today. This is natural in our model where baryons, CDM 
and photons are all generated by both inflaton and curvaton decay In terms of the 
spectrum (13 7|) . this means f$ ~ 1, which is the generic value. To get a much larger or 
smaller value, one has to tune the model in some way (say, by having A 7 = A C( i m = Ab 
to a high accuracy). 

Uncorrelated adiabatic and isocurvature perturbations add power on the SW 
plateau relative to the acoustic peaks, and since the observations indicate a deficit of 
power, this does not improve the fit to the data jH] . However, in our model the adiabatic 
and isocurvature perturbations are correlated. Negative correlation^ will naturally 
suppress the power on the SW plateau and thus improve the fit. The suppression 
is large if TZ ~ S, and total when 1Z = 2(1 — f v )S. Again, note that our model 
naturally yields TZ ~ S which is required for significant cancellation. Also, as we will 
discuss in section our spectrum (jHTjl gives more freedom to adjust the cancellation 
of the low multipoles than in previous studies of suppression due to anticorrelation 

[ainiiniEiiEiiEHiEnniniEiiioi]. 

We note in passing that the conclusion that the peak structure excludes a sizeable 
contribution from CDM isocurvature modes outside of the SW plateau does not hold for 
neutrino density and neutrino velocity isocurvature modes, since their peak structure is 
similar to that of adiabatic modes [1011221 • In fact, there is an explicit example where 
isocurvature modes give a major contribution to all of the observed acoustic peaks 

3.3. Fitting method 

We fit the spectra ® and flU to the WMAP first-year data ^ EH E3 by running 
several Monte Carlo chains on a supercomputer. We vary the following 11 primary 
parameters (parameter ranges are given in Table Q). We have five standard parameters 
related to the cosmological background, called the hard parameters: the physical baryon 
density ujf, = h 2 Qb, the physical dark matter density ujdm — h 2 Qdm, the ratio of the sound 
horizon to the angular diameter distance (i.e. the acoustic peak scale in radians) 0, 
the optical depth due to reionisation r and the neutrino fraction f v = u u /ud m . We also 
assume that there is vacuum energy, the density of which is fixed by the constraint that 

|| The factor 1 — /„ arises because we are considering the cold dark matter-photon isocurvature 
perturbation S c d m ,~f instead of the total dark matter-photon isocurvature perturbation Sdm~ r This 
issue will be discussed in section ET5l 

% Note that in our convention negative correlation for the Cy.s on the SW plateau corresponds to 
positive correlation for TZ and S, since g[ s {k) = —2(1 — fu)gj-ji{k). One should be careful, as sign 
conventions vary. Some authors prefer to use C = —TZ, and others define the variables with a different 
sign; e.g. in |44 | there was an additional minus sign in the definition of TZ, so that negative correlation 
between TZ and S corresponded to a negative C i cor . 
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the universe is spatially flat, {l\+£l m = 1, where Q m = fi&+f^ m . We have six parameters 
for the spectrum of perturbations, instead of the usual two: the spectral index m, 
the difference between the spectral indices An = ri\ — ri2, the overall normalisation 
ln(10 10 A^), the fraction of adiabatic perturbations that come from the inflaton cos(/>, 
the fraction of isocurvature perturbations that come from the inflaton cos#, and the 
CDM-photon isocurvature fraction f$- 

The technical details of our analysis are very similar to those described in the 
appendix of |HBj- We use a modified version of CosmoMC [HE] to generate 40 Monte Carlo 
Markov Chains that start from randomly chosen points in the parameter space. The 
transfer functions gf n (k), gf s (k), gfn(k), and gf s {k) are calculated with our modified 
version of CAMB jHO] whenever any hard parameter is changed. The integrals 
and (J40|) are then evaluated for each model. After running between 8 and 12 days, 
each chain contains about 10 000 accepted steps (150 000 step trials). We cut off 
burn-in periods and neglect 13 chains that never burned-in but instead got stuck in a 
local likelihood minimum. The correlation length (i.e. the number of steps between 
independent samples) is quite long in our chains, typically 170 - 1000. Nevertheless, 
after the described procedure, we still have over 220 000 accepted steps (3 300 000 step 
trials) giving 7500 independent samples to analyse. 

In the same way, we create chains for a ACDM model with adiabatic power-law 
perturbations (and without massive neutrinos) to serve as a reference point. The 6 
primary parameters of the reference model are u>b, Udmi ® an d t for the background and 
n, ln(10 10 A^) for the perturbations. 

We produce a 1-dimensional marginalised likelihood distribution (probability 
density) for each primary parameter and for some derived parameters such as the spectral 
index n 2 , the dark matter-photon isocurvature fraction f$, the vacuum energy density 
Q\, the baryon energy density Qf,, the matter energy density Q m , the CDM energy 
density fl c dm, the (massive) neutrino energy density Q u , the age of the universe, the 
reionisation redshift z re and the Hubble parameter today H = h x 100 km/s/Mpc. 

We do not apply any priors to primary parameters, but use a top-hat prior 
0.21 < h < 1.00 for the Hubble parameter and Q\ > for the vacuum energy density. 
The main focus of the present paper is the study of the perturbation spectrum in the 
mixed inflaton-curvaton model, so we will here simply do a fit to the WMAP data, and 
briefly discuss the interesting features of the fit. We will follow up with a more detailed 
analysis of cosmological parameter evaluation in a separate paper j2S], where we will 
also include data from high-/ CMB measurements and the Sloan Digital Sky Survey 

(sdss) HHEI]. 

3.4- Results of the data fit 

There are 1348 WMAP data points and our model has 11 parameters, so the reduced 
number of degrees of freedom is v — 1348 — 11 = 1337. The reference pure adiabatic 
ACDM model has 6 parameters giving v = 1342. Our best-fit model has x 2 = 1423.9 
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and the best-fit reference model has x 2 = 1429.0. In the first two columns of Table 
[0 we compare the parameters of these models. Even though our model has a better 
X 2 value, x 2 / u turns out to be exactly the same (1.065) for both models. Thus this 
measure-of-goodness does not support introducing the isocurvature degrees of freedom" 1 ". 
However, we are interested in the SW plateau (in particular the quadrupole power) 
beyond its statistical weight in the overall fit. On the last three lines of Table ^ we 
give the quadrupole power of both models as well as the x 2 °f the first two data points 
(quadrupole and octopole) and their contribution to the total x 2 ■ The quadrupole power 
in our best-fit model is 125 /iK 2 smaller than in the ACDM model. This makes our fit 
better in the SW plateau and drops the x 2 contribution of the first two data points from 
0.51% to 0.45%. 

The parameters of the best-fit model are somewhat exotic. The physical baryon 
density cj& = 0.041 is much larger than the value given by BBN |SH]. The large neutrino 
density f u = 0.86 implies (along with uj^m = 0.12) that the sum of neutrino masses 
is 9.2 eV, in conflict with measurements of tritium decay which give an upper limit of 
m Ue = 2.2 eV for the electron neutrino mass [HI], yielding an upper limit of 6.6 eV for the 
sum of neutrino masses (when combined with the small mass differences from oscillation 
experiments). Note also the very small CDM density, £l c d m — 0.028. One would expect 
this to be highly problematic for structure formation, but on the other hand the spectral 
index on the relevant scales is deeply blue, rii = 3.9, which significantly enhances the 
amplitude of perturbations on the scales probed by the measurements of large-scale 
structure. 

The best-fit isocurvature model (for which we do not give a Ci plot) is not meant 
to be taken as an alternative to the best-fit adiabatic power-law ACDM model. Rather, 
it underlines the point that models with radically different background cosmologies and 
perturbation spectra can fit the WMAP data even slightly better than the standard 
model. Including CDM isocurvature perturbations opens up the space of well-fitting 
models, and the cosmological parameters in these models are no longer clustered around 
the best-fitting model. The addition of priors from other cosmological observations, as 
well as data on the power spectrum from high-/ CMB measurements and large-scale 
structure will exclude many of these models. For example, our best-fit model does not fit 
the data from high-/ CMB measurements: the nearly scale-free spectral index ri2 = 0.988 
dominates until the third acoustic peak, but the large spectral index rii = 3.906 gives 
too much power in the region / > 900. 

For now, the message to take away is that one cannot determine the cosmological 
parameters from the WMAP data alone with any degree of confidence without making 
strong assumptions about the primordial power spectrum. Turning this around, it is 
not possible to determine the primordial power spectrum using only the WMAP data 
without making strong assumptions about the cosmological parameters. This has been 
previously emphasised in [SI3 ESI- As a warning example, we note that the original 

+ As discussed in |fi2j . the criteria for introducing new parameters should in fact be more strict than 
simply getting a slightly better x 2 /v. 
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implementation of the pre-big bang scenario [fJ7] was ruled out in part because the 
spectral index n = 4 of the adiabatic perturbations coming from the dilaton was too 
large to fit the CMB data, and the second version jHH| was ruled out because the axion 
with a spectral index n = 1 carried only isocurvature perturbations. Implementing the 
pre-big bang curvaton model for the axion 0^3' but keeping the dilaton perturbations 
would lead to our spectrum (JHTj) with the spectral indices rii = 4 and n 2 = 1, which 
incidentally agree with our best-fit model. 

3.5. Likelihood distributions 

Let us briefly discuss the main features of the marginalised likelihood distributions, 
concentrating on the differences between our isocurvature model and the pure adiabatic 
ACDM model. 

In Table ^ we give the median values for the primary and derived parameters, and 
the regions that contain 68% of the accepted models in the 1-dimensional marginalised 
likelihood distributions for the adiabatic model and for our model. As the likelihood 
function is non-Gaussian with respect to some parameters, most of the quoted numbers 
cannot be interpreted as la confidence levels and are for comparison purposes only. (The 
meaning of the given values is as described in the appendix of [EE!-) Of the primary and 
derived parameters, only u^m, @, ^2, ^cdm and the age of the universe have a strongly 
Gaussian 1-d marginalised likelihood distribution. Of course, even for these parameters, 
the numbers are confidence ranges, and should not be treated as exclusion limits, as 
emphasised in jHE]- Adding more input into our analysis, such as the SDSS data [qTHIBT] 
would lead to better behaved likelihood distributions, in agreement with the discussion 
above. 

The well-known degeneracy between the physical baryon density and 
isocurvature modes |4^| EH] is apparent also from our results. In a pure adiabatic 
power-law model, is determined simply by the relative heights of the successive 
acoustic peaks. Since the CDM isocurvature mode has valleys where the adiabatic 
mode has peaks, increasing the relative isocurvature contribution can compensate for 
adding baryons. This, and the fact that the modes are the sum of two power-laws, 
means that there is much more freedom for the values of uj^ in our model than in a pure 
adiabatic model. Also, the preferred values of oj^ with isocurvature modes are generally 
higher than in the pure adiabatic case. 

The optical depth r is poorly constrained in the pure adiabatic case, and even more 
so in our model. The 95% region spans a huge range from 0.05 to 0.59, and the 68% 
region given in Table ^ is also wide. The optical depth is primarily determined by the 
single data point at I — 2 in the WMAP TE data. This quadrupole TE power is very 
high, hinting at early reionisation (large optical depth) in the pure adiabatic case. In 
our model the mechanism that suppresses the low TT multipoles also cancels TE power, 
which means that an even larger r is needed to fit the high first data point; on the other 
hand, positive correlation would add TE power at the quadrupole and allow a smaller 
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4.031+°^ 


cos 4> [0, 1] 


(1) 


0.295 


0.422 


(1) 


no constraint 


cos<9 [-1, 1] 




0.985 


0.731 




> 0.76 


fs [-42, 42] 


(0) 


-9.156 


-1.400 


(0) 


\fs\ < 13 


Derived parameters 












^2 




0.988 


0.845 




0.869+n 1 - 2 ^ 

— U.lou 


/5 = (1 - /,)/5 


(0) 


-1.329 


-0.385 


(0) 


\fs\ < 0.7 


n b 


0.695 
0.048 
0.305 


0.736 
0.069 
0.264 


0.345 
0.103 
0.655 


+0.087 
u -' ZD -0.084 

0.046+^08 

u 1 —0.087 


627+ - 141 

U - DZ ' -0.243 

0.079 + H?q 
0.373+°- 24 f 

u.u i ^_o.l41 




0.257 


0.028 


0.151 


0.229+g-gS 
— u.u / y 


< 0.06 


0^ 




0.166 


0.399 




0.239 + n?o4 


Age (Gyr) 

z re 

H (km/s/Mpc) 


13.62 
13.73 
68.87 


12.68 
25.93 
77.02 


14.32 
13.45 
52.21 


13.39±g;|S 
16.861^ 

71 70+lOdO 


13.691?-^ 

i8-4i+?:? 8 7 

D *- zy -11.38 


Low-I behaviour 












2(2 + l)C 2 /2vr (fiK 2 ) 


1170 


1045 


736 


(mean) 1262t 2 9 l : 


(mean) 1081^ 


x 2 (c 2 ,c 3 ) 


7.338 


6.408 


4.336 


(mean) 7.888 


(mean) 6.615 


x 2 (c 2 ,c 3 )/ x 2 (%) 


0.514 


0.450 


0.304 


(mean) 0.550 


(mean) 0.461 



Table 1. Comparison of models. In the last two columns we give the median of 
the 1-dimensional marginalised likelihood for each parameter and the region around 
the median that contains 68% of the accepted models in our chains. If the likelihood 
function was exactly Gaussian, this would correspond to the la region. However, most 
of the likelihood functions are too far from Gaussian for this to be true, as should be 
clear from the values of our best-fit model. Thus the parameter ranges should not be 
taken as ler regions. They are given for comparison purposes only. 

t. Note that in the pure adiabatic case, the signal for the large optical depth in the 
TT spectrum comes completely from the southern hemisphere; the preferred value for 
the northern hemisphere is r = Such asymmetry suggests caution about the 

interpretation of the data points leading to the value of r. At any rate, the effect on 
our results is small, since they are consistent with a wide range of values. 

In the pure adiabatic power-law model there is only one spectral index while our 
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model has two, so a simple comparison of their values is not very meaningful. However, 
let us mention that while the peak of the likelihood function for the larger spectral index 
ni is at the moderate value 1.1, the distribution has a long tail at large values, so that 
the median is 1.4 and the mean is even larger, 1.8. While this may be an artifact of poor 
convergence of An (the split test shows that the upper bound for An is not reliable), 
our best-fit model with rti = 3.9 demonstrates that the WMAP data does not require 
the spectrum to be almost scale-invariant over the entire range of scales covered by the 
data. 

The vacuum energy density is much less constrained than in pure adiabatic 
models. The reason is again low-/ modifications. Increasing Q\ raises the TT power 
at the lowest multipoles due to the integrated Sachs- Wolfe effect. This effect can be 
countered by negative correlation which brings the low multipoles down, or mimicked 
by positive correlation which takes them up. The net effect is to widen the range of 
Qa in both directions. The Hubble parameter and vacuum energy density are strongly 
correlated, so the freedom in the values of is reflected in the values of Ho- 
lt is notable that our isocurvature spectrum prefers a large neutrino fraction /„, 
in contrast to a pure adiabatic spectrum which is not very sensitive to the value of 
f v [HE EDI EH 172] • Our model differentiates between CDM and massive neutrinos 
because the latter carry no isocurvature perturbations. This means that for fixed dark 
matter density and CDM isocurvature amplitude, one can tune down the isocurvature 
contribution by increasing the neutrino fraction f v : the less CDM there is, the 
less difference its perturbations make. The total dark matter-photon isocurvature 
perturbation Sd m ^ and cold dark matter-photon isocurvature perturbation iS c d m)7 are 
related by Sd m ,-y = (1 — f^Scdm^, so the power spectrum is almost (but not entirely) 
degenerate with respect to scaling f s and 1/(1 — f u ) by the same factor. As it is, our 
model seems to prefer a high CDM-photon isocurvature fraction fs, with a large neutrino 
fraction to compensate. A meaningful measure of the isocurvature contribution is given 
by the dark matter-photon isocurvature fraction fs = (1 — fu)fs- We find \ f$\ < 0.7 
in the 68% region (in agreement with j3H EHl EH] ) , though a look at our best-fit model 
with fs = —1.3 should again serve as a warning not to interpret this number as an 
exclusion limit. 

Since fs and 1 — f v are almost degenerate, measuring one will constrain the other. 
As CDM and neutrinos have a very different effect on the matter power spectrum [58] l7Uj. 
using the SDSS data to give a handle on the neutrino fraction could indirectly constrain 
the isocurvature contribution (which is itself negligible at the SDSS scales since the 
isocurvature perturbations are damped rapidly with increasing /). 

Finally, from Tabled we see that the mean TT quadrupole power in the MC chains 
for our model is 181 uK 2 smaller than for pure adiabatic models. This shows that 
the cancellation mechanism works and on average helps in fitting the low-/ part of the 
spectrum. The mean x 2 contribution of the first two TT data points to the total x 2 is 
also reduced from 0.55% to 0.46%. We devote the next section to this important issue. 
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3. 6. Suppression of the low multipoles 

The most important feature of the angular power spectrum that distinguishes the 
inflaton-curvaton model from the standard adiabatic power-law ACDM model is the 
amplitude of the low multipoles. In fact, this is the only significant qualitative difference, 
which is not surprising. As discussed in section 13 .2^ the pure adiabatic ACDM model 
fits the TT power spectrum very well, apart from the low multipoles and the three 
"blips" at intermediate multipoles before and at the first acoustic peak (see Figure EJ). 
The three blips cannot be fitted with the smooth CDM isocurvature modes, so the only 
place where we can get significant improvement is the low multipoles. 

While the discrepancy of the low multipoles in the adiabatic ACDM model looks 
significant to the eye (see Figure EJ), the small number of badly fitted data points 
(practically the first two points) combined with the large cosmic variance means that 
the weight of this feature in the overall x 2 fit is quite small (for the best-fit adiabatic 
model only 0.51%, or 7 points, of the total y 2 ). Therefore one does not expect any 
model that does not explain the "blips" to have a significantly lower y 2 than the pure 
adiabatic ACDM case, regardless of what the model is. Indeed, the difference in x 2 
between the best-fit pure adiabatic model and the best-fit inflaton-curvaton model is 
only 5: our model fits better, but not decisively better. 

However, the small amplitude of the low multipoles is an interesting feature beyond 
its statistical significance, and has attracted a lot of attention. Since a x 2 -selection 
among models does not particularly pick out ones with small amplitude for the low 
multipoles (in our best-fit model, less than 1 point of the improvement in the x 2 comes 
from the quadrupole and octopole), we specifically look for low quadrupole power among 
those of our models that fit the data nearly as well as the best-fit model. (The inclusion 
of data from high-Z CMB measurements and SDSS would further reduce the weight of 
the low multipoles in the overall fit, making it more difficult to find models with a low 
quadrupole.) 

We take a model found using this method as our example with which to demonstrate 
the suppression mechanism of the low multipoles. The parameters for the example model 
are given in the third column of Table ^ The TT and TE angular power spectra of the 
example model and the contributions from the individual adiabatic, isocurvature and 
correlation components are shown in Figure El In Figure E] the inflaton and curvaton 
contributions and the total TT spectrum are plotted in comparison to the adiabatic 
power-law ACDM model. 

The example model has x 2 = 1425.9, which is 2 larger than our best-fit, and still 3 
smaller than the adiabatic best-fit. On the last three lines of Tabled we compare the 
best-fit adiabatic ACDM model, our best-fit model and our example model. We give 
the quadrupole power, the x 2 °f the first two data points (quadrupole and octopole) 
and their contribution to the total x 2 ■ The quadrupole power in our example model is 
736 /uK 2 , which is 434 /iK 2 smaller than in the best-fit adiabatic model. This makes 
the fit better on the SW plateau and drops the x 2 contribution of the first two data 
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points from 0.51% to 0.30%. Note that unlike in our best-fit model, practically all of 
the improvement in the x 2 over the adiabatic case comes from the low multipoles. 

For comparison purposes, we also did a search for low quadrupole power among 
ACDM models that are within A% 2 < 2 from the best-fit ACDM model. As expected, 
the lowest quadrupole power found (914 //K 2 ) is significantly higher than in our example 
model. 



Extreme cancellation. Let us look at the details of the suppression mechanism of the 
low multipoles with the spectrum (J37j) . If we wanted to pull the amplitude on the SW 
plateau down as much as possible, we could put cos0 = in 



to kill Cf dil and 



Cf orl , and then use anticorrelation to get rid of (7f dl2 + Cf° 2 + Cf or2 on ^ e SW plateau. 



Given that gf s {k) = —2(1 — fv)gJn{k) on the SW plateau, this implies, from (|HHJ), that 
sin</> + 2(1 — f„)fssm6 = 0. The TT spectrum would then be 

cot 2 ^"- 1 + ( Q Uk) + 2 9 ^ k) f) \ ~k n ^ 



TT 



Air 



— N 

k 



K 



2(1 -/„) 



(41) 



and we have cos</> — 0,fs — — 1/(2 sin 6). On the SW plateau, the second term is zero, 
and the spectrum has the index ni, with an amplitude that is freely adjustable with 
cot 2 8. Away from the plateau, the isocurvature transfer function gf s {k) is damped more 
rapidly than the adiabatic transfer function gf-ji(k), so the spectrum at high multipoles is 
essentially adiabatic, has the spectral index n 2 and an amplitude given by N\. Since we 
have ni > n 2 , the low multipoles are further suppressed relative to the high multipoles. 

In other words, the low multipoles are given by inflaton isocurvature perturbations 
with the larger spectral index ni, while the peak structure is given by adiabatic curvaton 
perturbations with the smaller spectral index n 2 . The SW plateau and the peak 
structure can therefore be adjusted almost independently to fit the data. 



Cancellation in our example model. In contrast to the extreme case, in our example 
model the behaviour of the high multipoles is not completely decoupled from that of 
the low multipoles. However, the cancellation mechanism for the low multipoles works 
qualitatively as sketched above: a small cos</> suppresses the adiabatic and correlation 
components C ; adl1 and C z corl with the larger spectral index ri\, and strong anticorrelation 
suppresses all three components Cf dl2 , Q so2 and C ; cor2 with the smaller spectral index 
n 2 . 

From FiguresElandiniwe see that the curvaton perturbations with the spectral index 
n 2 = 0.85 are killed almost completely by the anticorrelation at the low multipoles, but 
they dominate the peak structure. The power at low multipoles comes almost completely 
from the inflaton perturbations, with ri\ = 1.07. However, in contrast to the extreme 
case, the amplitude of the adiabatic component Cf dl1 is not exactly zero, and it gives 
a significant subdominant contribution to the acoustic peaks. For the first peak, the 
inflaton and curvaton correlation modes C, corl and C, cor2 also contribute. The correlation 
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Figure 5. The example model of the third column of TablcHwith the WMAP data (•) 
and the highest I data points of other CMB experiments from Tcgmark's compilation 
[7D| (♦). The total angular power consists of the six components mentioned in the 
paragraph after equation (gO): Cf ot = Cf dil + Cf ol + C*f orl + Cf di2 + Cf ° 2 + C*f or2 . 
Upper panel: Temperature-temperature angular power 1(1 + l)Cf T /2tt. Lower panel: 



Temperature-polarisation cross-correlation power (I 



l)Cf E /2K. 
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Figure 6. Comparison of our example model (the third column of Tabled see a ls° 
Figurc|SJ) to the best-fit 6-parameter adiabatic ACDM model (the first column of Table 
GJ. The total angular power of our model is the black solid line; C* ot 
where C*f uml = Cf dil + Cf ol + Cf ml and Cf um2 = Cf di2 
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results in peak structure between pure adiabatic and pure isocurvature cases, and the 
fact that those components are themselves a sum of two power-laws also slightly modifies 
the relative heights of the successive peaks (and valleys). 

Since our example model has a non-zero vacuum energy density (Q\ = 0.35), 
the integrated SW effect would typically lead to rising power towards the smallest 
multipoles (as apparent for the adiabatic ACDM model in Figure EJ). However, the 
inflaton perturbations are going down with increasing / because they are dominated 
by the isocurvature component while the curvaton perturbations are slowly rising with 
increasing I as the cancellation by anticorrelation weakens after I ~ 10 (as the SW 
effect ceases to be dominant), and this interplay keeps the spectrum quite flat at low 
multipoles. 

Not all features of our example model are universal to well fitting models with 
the spectrum (|3*7f). As discussed in section EP1 good models do not cluster around a 
single region in parameter space. However, the suppression mechanism which combines 
anticorrelation and two different spectral indices is indeed a generic feature of our 
spectrum (|37[). and is not tied to the exotic features (such as the high physical baryon 
density of our example model. 
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Running of the spectral index. Analysis of CMB data with pure adiabatic models with 
a running spectral index has pointed towards the feature that the spectral index evolves 
from a large to a small value as k increases |HI1 EH ESj , though the evidence is marginal. 
(Note also [7^] where the fit to the data is improved by having an index that does not 
run but changes discontinuously from a large to a small value.) 

For a pure adiabatic spectrum which is a sum of two different power-laws, the 
effective running of the spectral index is always positive: the smaller index dominates 
at small k and the larger index at large k. As discussed above, in our model it is natural 
for the large index to dominate at small k and vice versa. Of course, suppression of the 
low multipoles via anticorrelation itself also looks like negative running of the spectral 
index if one tries to fit it with a pure adiabatic power-law. Indeed, the indications of 
running in the WMAP data come mostly from the amplitude of the low multipoles |74j . 

Correlated isocurvature perturbations allow negative running because the prefactors 
evolve with k, as seen in (J39|) . At small k, the prefactor of the perturbations with 
the small index is suppressed by anticorrelation, while at large k the prefactor of the 
perturbations with the large index is suppressed by the rapid damping of isocurvature 
perturbations with k. A key feature allowing the complete cancellation of the adiabatic 
terms at the low multipoles and thus the decoupling of the low and high multipoles is 
that the adiabatic and isocurvature spectra have the same shape, a combination of two 
power-laws. 

Previous studies where the adiabatic and isocurvature spectra have been taken to 
obey a simple power-law with the same index [T3 ^ l4T | l5U] have found similar suppression. 
However, that form of the spectrum is too constrained for the low and high multipoles to 
be decoupled. For more complicated spectra, the adiabatic and isocurvature components 
have often not had the same shape EH EHl EH E2] , so they cannot effectively cancel, 
and the low multipoles again cannot be as decoupled from the high ones and as effectively 
suppressed as in our case. (However, double inflation models with spectra similar to 
(|37|) have been studied |2~T) I22|.) In the case of models with more than one isocurvature 
component, but all having the same spectral index [SH ESI ESj , the fit to the data has 
also improved less over the adiabatic case than for our spectrum (|37|). Note that this is 
not a question of the number of parameters: for example, the models studied in [J2JEI] 
have as many or more parameters than our model. This underlines the importance of 
priors in the form of the spectrum on the estimation of the isocurvature contribution. 

4. Conclusion 

We have studied a model that generalises the usual curvaton scenario in two ways. First, 
we model the curvaton properly as a scalar field instead of dust when it decays. We 
follow the curvaton decay numerically. We find that the results of the dust case are 
recovered whenever the curvaton does not significantly contribute to the energy density 
when it is starting to oscillate, even if the curvaton does not oscillate rapidly when it 
decays. 
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Second, we take into account that in addition to the curvaton, the inflaton can 
also have perturbations. The resulting spectrum is an interesting mixture of correlated 
adiabatic and CDM isocurvature perturbations, which could arise also from a double 
inflation model. We analyse this spectrum using the first year WMAP data. We find that 
including the isocurvature modes opens up the parameter space of well fitting models, 
and the parameters for both the perturbations and the cosmological background no 
longer cluster around the best-fit of the standard adiabatic power-law ACDM model. 
For example, the spectral indices do not need to be close to scale-invariant, and the 
physical baryon density ui, is typically much larger than the usual value. Adding other 
input, such as data from high-/ CMB measurements and SDSS would exclude many of 
the new well-fitting models. We will consider these issues in a more thorough analysis 
of the perturbation spectrum in a follow-up paper [23] . 

The most significant qualitative feature of our model is the suppression of the low 
multipoles due to anticorrelation between adiabatic and isocurvature perturbations. Our 
spectrum makes it possible to decouple the behaviour of the low and high multipoles and 
thus fit the observed small amplitude of the quadrupole and octopole without affecting 
the adiabatic peak structure, improving the fit to the data. There is no need to fine tune 
the multipoles affected by the suppression: if the perturbations in the decay products 
of two (or more) scalar fields both contribute to the radiation and matter observed 
today, it is the low multipoles that are generically affected by the resulting isocurvature 
contribution. Recently, another model where the suppressed multipoles do not have to 
be picked by hand, but are instead related to the dark energy via an IR/UV duality, 
has been proposed [77j . 

The significance of having qualitatively different behaviour of the spectrum at small 
and large k goes beyond explaining the small amplitude of the low multipoles, because 
the determination of cosmological parameters is sensitive to the behaviour at small k. 
For example, in [ZEj a model with pure adiabatic perturbations with a different spectral 
index for small and large k and no cosmological constant was shown to fit the CMB and 
large scale structure data better than the ACDM model. 

There are indications that the WMAP signal for the quadrupole as well as 
other low multipoles, up to / = 40, has significant non-cosmological contamination 
[OHl EH1 EE EDI IHU E2j This could completely change the estimation of isocurvature 
perturbations, since their contribution is determined mostly from the low multipole 
part of the spectrum. On the other hand, the presence of such contamination would 
underline the importance of studying how cosmological parameter estimation is affected 
by possibly obscured features at small k, such as isocurvature modes. 
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